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ABSTRACT 

The motivation for Computational Aeroelasticity (CA) 
and the elements of one type of the analysis or 
simulation process are briefly reviewed. The need for 
streamlining and improving the overall process to 
reduce elapsed time and improve overall accuracy is 
discussed Further effort is needed to establish the 
credibility of the methodology, obtain experience, and 
to incorporate the experience base to simplify the 
method for future use. Experience with the application 
of a variety of Computational Aeroelasticity programs 
is summarized for the transonic flutter of two wings, 
the AGARD 445.6 wing and a typical business jet 
wing. There is a compelling need for a broad range of 
additional flutter test cases for further comparisons. 
Some existing data sets that may offer CA challenges 
are presented. 

INTRODUCTION 

One troublesome area for aeroelastic analysis has been 
the transonic speed range. Shock waves that occur in 
the flow over wings and bodies are not included in 
conventional linear theories. Furthermore, minimum 
flutter speeds, buffeting, limit-cycle-oscillations 
(LCO), aileron buzz, and shock-boundary layer 
oscillations may be encountered. Computational Fluid 
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Dynamic (CFD) methods have been extensively 
developed and applied in the area of steady aerodynamics 
in the past 2-3 decades. Application of this 
methodology for unsteady aerodynamics to be used in 
aeroelastic analyses has also been under development. 
When CFD is coupled with the structural dynamics in 
the computational process, it is generally referred to as 
Computational Aeroelasticity (CA). One of the 
difficulties with CA is that aircraft must be cleared for 
flutter significantly beyond cruise conditions where 
strong shocks, separated boundary layers, and other 
flow conditions which are difficult to compute, may be 
encountered. 

In the past, CA has required long solution times on 
expensive supercomputers. Thus, the associated costs 
were very large. However, in the past decade, 
workstation-type machines have attained the 
performance level of the supercomputers of the 
previous decade and the cost of the computation has 
decreased by between two and three orders of magnitude. 
Nevertheless, the cost of labor, elapsed time, and 
peripheral processing is still significant. The 
technology has evolved to the status that it can be used 
for analyses to check cases that might be considered 
critical based on linear analysis. However for 
multidisciplinary design problems which can involve as 
many as 10 5 to 10 6 evaluations of the flutter eigenvalue 
matrix for stability, CA is impractical without further 
development of innovative methods of application. 
This paper will review one approach to an aeroelastic 
analysis process with some comments on possible 
improvements. Some CA results will be discussed for 
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Fig. 1 Computational Aeroelasticity analysis block diagram 
for time -domain analysis. 


two wings, the AGARD 445.6 wing and a typical 
business jet. The oscillatory shock-boundary layer 
oscillation phenomenon will also be highlighted and 
some possible challenges for the CA researchers to 
consider are presented. 

OVERVIEW OF ONE CA ANALYSIS PROCESS 

One approach to computational aeroelastic analysis is 
depicted by the block diagram of Fig. 1. This 
illustrates a time-domain approach based on 
representing the structural dynamics with vibration 
modes. As indicated in the Introduction, emphasis 
should be placed on making the entire process more 
efficient. In this section some overall comments are 
given on the CA procedures in use at the Aeroelasticity 
Branch of the NASA Langley Research Center (LaRC). 
Similar approaches at other organizations are presented 
in refs. 1-3. 

In the Aeroelasticity Branch there are several CFD 
codes in use. At the transonic small disturbance level, 

4 

the in-house developed code CAP-TSD 
(Computational Aeroelasticity Program - Transonic 
Small Disturbance) is used. Also viscous effects can 


be considered with a stripwise, inverse boundary layer 
code incorporated in CAP-TSD by Edwards 5 and called 
CAP-TSDV. Codes at the Euler and Navier-Stokes 
level are the CFL3D code developed by the 
Aerodynamic and Acoustic Methods Branch at LaRC, 

7 

and ENS3DAE which was developed at Lockheed- 
Martin under Air Force sponsorship. 

The three major elements of most CA methods involve 
preprocessing of the geometry and modal data, the 
execution of the CFD computer program, and the 
postprocessing of the output data. Of course at all 
stages of the process, plotting and monitoring of all 
data is essential both for assuring that the data are error 
free, and that the results are properly converged. The 
plotting and monitoring tasks can be very time 
consuming and a high level of automation is desired. 

For mode processing which requires interpolation from 
the vibration-modal grid to the aerodynamic grid, we 

g 

have generally used the surface spline. For the small 
disturbance code that requires modal slopes as well as 
amplitudes, some of the limitations of this spline can 
be overcome for cases using calculated vibration modes 
by splitting the modal rotations rather than using the 
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slope of the spline. This capability is basically 
developed for wings only. The development of a 
suitable method for including fuselage or bodies may 
require significant effort. An evaluation of this 
technology is given in ref. 10. 

Geometry processing is usually treated on a case-by- 
case basis, and can range from a simple spline for 
airfoils to a more elegant Computer-Aided-Design 
model (CAD). The TSD methodology also requires 
slopes, which demands particular care. The significant 
task is in grid generation for the CFD portion of the 
process. With TSD codes, the task can be built into a 
short interactive program making use of experience. 
For the higher equation level codes the grid generation 
is more complex and requires a significant level of user 
skill. 

For the TSD codes, the surface velocity boundary 
condition is applied at a mean plane and thus the grid 
does not have to be moved with the deforming surface. 
However, at the Euler and Navier-Stokes level, the grid 
must move with the surface. Normally this is done 
with some type of deforming mesh algorithm. There 
are several types in use, such as a spring analogy, 

7 

simple shearing, or more complex deforming 
algorithms exemplified by the complex multi-block 
case of ref. 12. Deforming meshes can be difficult as 
in certain cases the mesh can fold over on itself creating 
a negative computational cell volume and program 
failure. Also some algorithms require a significant 
increase in memory and CPU time over the basic CFD 
code. This is an active research topic at the present 
time. 

For this type of problem, CPU speed, algorithm 
efficiency , and parallel processing to decrease turn- 
around time are essential. The rule of thumb developed 
in some of the early CA flutter analyses is that it takes 
about 100 runs to calculate a flutter boundary for 
several Mach numbers and dynamic pressures using 
time simulation. This would include runs for static 
initial conditions, debugging, and dynamic response. 

For post processing of output time domain signals for 
linearized stability, the damped sine wave fitting 
procedure of ref. 1 3 has been extensively used. It is 
limited to analyzing a single modal coordinate and to a 
low number of modes in the fit, usually about six. 
Other methods have been developed along the lines of 
ref 14. In addition, some initial trials of the system 
identification of ref. 15, which is incorporated in the 
commercial program of ref. 16, appear quite 


promising. For low amplitude stability assessment, 
the record length requirements of the system 
identification method dictates the length of the 
computer run and is a strong driver on computer time. 

It might be noted that for symmetric motions only a 
half airplane needs to be treated by CA methods. 
However for antisymmetric motions, the full airplane 

must be treated thus doubling the memory and CPU 
17 

requirements. Furthermore, both symmetric and 
antisymmetric modes must be retained in the analysis 
because the static aeroelastic deformation for the initial 
condition involves the symmetric modes. The cost of 
the dynamic analysis is thus doubled as well unless 
there is special treatment of the deformed shape. 

The CA methodology requires both a reasonable 
appreciation of the CFD technology and of 
aeroelasticity. Generally it has been found to work best 
with a teaming approach. 

VERIFICATION 

Unfortunately it is relatively easy even for an 
experienced user to run CA programs and generate 
reasonable looking results that may not be accurate. 
This puts a burden on the user and developer to verify 
both the program and the results. A recent issue of the 

AIAA Journal contains a 100-page section on credible 
18-29 

CFD calculations " with emphasis on separating the 
issues such as grid convergence, residual convergence, 
algorithm convergence, and program correctness. 
Unfortunately, in the quest for time accuracy 
verification, there are very few test cases other than 
one-dimensional shock problems and one-dimensional 
wave propagation problems. With a TSD code it is 
relatively easy to revert back to linear theory for 
comparison, whereas gridding issues make it more 
complex for the higher level codes. Generally it is 
preferable to start a flutter analysis by verifying a result 
obtained from linear theory. 

Comparisons with experiment are essential for 
verification, but agreement with experiment does not in 

itself serve as complete program verification. 
30 

Bobbitt has discussed the CFD validation process and 
listed ten items in each of the experimental and CFD 
areas that may affect the outcome of comparisons. The 
list was primarily oriented toward static results. For 
dynamic or flutter cases, modal accuracy and 
completeness becomes a further factor. In addition 
flutter conditions measured in a wind tunnel are usually 
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are a judgment call by the investigator and should not 
be interpreted in the same sense usually inferred for 
static wind tunnel measurements. Striving for exact 
agreement between theory and experiment is futile. 
Sometimes there can be systematic errors such as 
indicated by one case of a comparative tunnel tests for 
unsteady data that surfaced an unusual effect of a root 
sidewall porosity that might not have been detected 
without the other tests.' 


AEROELASTIC TEST CASES 

As previously mentioned, an essential part of the 
verification process is the availability of good 
experimental data for test cases. One flutter test case 

that is publicly available with enough information for 

32 

CA application is the 445.6 wing which is known as 
AGARD Configuration I for which our results will be 
subsequently discussed. This configuration has been 
extensively used for computational studies. A limited 
survey yielded 15-20 different investigations using this 
data. Although it has been over ten years since this 
data set was made available, there is still no 
Configuration II. Many of the flutter data sets do not 
have adequate information for geometry or adequate 
modal information for benchmark calculations. Some 

efforts have been made in the LaRC Benchmark Models 

33 

Program where the data are for simple planforms 
only. During these tests, significant effoit was placed 
on measuring unsteady pressures during flutter. Some 
of the results will be discussed in a later section. 

Another approach is to try to calibrate the unsteady 
aerodynamics separately, by measuring unsteady 
pressures and/or forces during forced oscillations such 
as pitch, plunge, or control surface oscillations. There 
is a wide range of data available for this purpose. One 
notable effort by AGARD resulted in an organized set 
of test cases for several two and three-dimensional 
configurations and are published in ref. 34-35. These 
cases have been widely used for evaluating unsteady 
aerodynamics. There is currently an effort to assemble 
another set of cases and another document by a 
Working Group under the Applied Vehicle Technology 
Panel of Research and Technology Organization 
(formerly AGARD). These cases will also include 
cavity flows, dynamic stall, vortex flows, as well as 
flutter-oriented data sets. Publication is expected in the 
latter part of 1999. 

One of the difficulties with just looking at unsteady 
aerodynamic comparisons is that it is difficult to 


interpret the effect of imperfect agreement of pressures 
on aeroelastic problems. That is, how much change in 
flutter velocity would be produced by a deviation in 
pressure at a point on the wing? For example, one 
would normally expect more aeroelastic effects from 
deviations at the tip than otherwise, but each case must 
be treated individually. Technically this is a flutter 
sensitivity problem which is straightforward for a 
specific case, but it is not generally investigated. 

There is a general need for flutter test cases for a variety 
of configurations for further validation efforts. Further 
efforts in collecting and organizing existing data sets is 
needed, and further tests for this purpose need to carried 
out. 


REVIEW OF SOME RECENT RESULTS 

The results for two flutter test cases are presented. The 
results are mostly those generated by the Aeroelasticity 
Branch and also have been presented in ref. 36-37. 
The flutter stability boundaries are presented in terms of 
flutter speed index, which is the flutter velocity divided 
by the root semichord and the square root of the mass 
ratio, and the frequency ratio which is the flutter 
frequency divided by the torsional frequency. 

445.6 Wings 

The 445.6 wing planform is shown in Fig. 2. There 
were several of these semispan models which were 
cantilevered from the wind tunnel side wall and had a 



Fig. 2 Planview of AGARD Wing 445.6 
Standard Aeroelastic Configuration 
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quarter-chord sweep angle of 45 deg (leading edge sweep 
of 46.3 deg), a panel aspect ratio of 1.65, and a taper 
ratio of 0.66. This series of wings was flutter tested 
both in air and in heavy gas in the Transonic Dynamics 
Tunnel (TDT) at the NASA Langley Research Center. 

q experiment 

CAP-TSO. Ilneai 9 



(a) Flutter speed index. 



Mach Number 

(b) Frequency ratio. 

Fig. 3 Comparison between experimental 
and calculated flutter speed index and 
frequency for the AGARD Wing 445.6 tested 
in air. 


The wings had a NACA 65A004 airfoil and were 
constructed of laminated mahogany. To reduce the 
stiffness, some wings had holes drilled through them 
and filled with foam. The vibration modes were 
calculated and are published in ref. 32. 

Fig. 3 summarizes the flutter results for the 445.6 
wing tested in air in terms of the flutter speed index and 
frequency ratio. For these cases, the mass ratio is 
relatively high, 0(100), and the reduced frequencies in 
the transonic range is relatively low, that is, less than 
0.10. The linear theory results from CAP-TSD, as 
indicated by the dashed line in Fig. 3, agree quite well 
with the measured values except in the range of Mach 
numbers near 1 .0 where the results are somewhat high. 
For this thin wing, transonic effects begin quite close 
to a Mach number of 1.0. Including thickness in the 
CAP-TSD results lowers the flutter speed slightly for 
the subsonic Mach numbers with the most reduction at 
M=0.96. The addition of the boundary layer in CAP- 
TSDV slightly raises the flutter speeds and they are in 
agreement with the experimental data. The CFL3D- 
Euler ' results are in good agreement at the lower two 

Mach numbers, but are increasingly low near M=0.96. 

39 

The CFL3D-Navier-Stokes result at M=0.96 is 
somewhat higher and shows improved agreement. 

For supersonic speeds, the CFL3D-Euler results show a 
premature rise and a large overprediction of the flutter 
boundary. Similar results have been shown in ref 40, 
and unpublished results with CAP-TSD, but the Euler 
results of ref. 41 are close to the experimental data. 
Here, the one point at M=1.14 with CFL3D-Navier 
Stokes shows a large reduction of flutter speed with the 
inclusion of the boundary layer effects, but is still 
about 18 percent high in V (or 40 percent in dynamic 
pressure). Further effoit is needed to investigate this 
sensitivity to viscous effects. 

Fig. 4 presents the CAP-TSDV results for the weak 
wings in heavy gas. The mass ratios for these cases 
run from 12-34 and reduced frequency from 0.33 to 0.18 
as Mach number is increased. Generally good 
agreement is shown in this case also. A dip calculated 
with CAP-TSDV is shown which occurs just before 
the last measured point near M=1 .0. Reynolds number 
and amplitude effects were explored showing a 

small Reynolds number effect near the dip and a 
moderate effect of excitation amplitude on the boundary 
at Mach numbers higher than the dip. 
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Mach number 

(a) Flutter speed index. 



Mach number 


(b) Frequency ratio. 

Fig. 4 Comparison between experimental 
and calculated flutter speed index and 
frequency for the AGARD Wing 445.6 tested 
in heavy gas. 

Business Jet Wing 

The business jet flutter model is shown in Fig. 5 
mounted in the Transonic Dynamics Tunnel. An 
extensive set of flutter and static aeroelastic calculations 
using CAP-TSD and CFL3D-Euler and Navier-Stokes 


were made by Gibbons ” to compare with the measured 
flutter boundary. Calculations using CAP-TSDV have 
been presented by Edwards. 



Fig. 5 Business jet flutter model mounted in 
NASA Langley Transonic Dynamics Tunnel. 

The model was 4.4 feet in semispan and was 
constructed from an aluminum plate with foam and 
fiberglass used to provide the airfoil shape. The 
flexible wing and the stiff fuselage were mounted to the 
turntable in the wall of the tunnel. The wing root 
angle of attack was varied to minimize loading. The 
maximum angle used for this purpose was 0.2 deg at 
the highest Mach number and this value was used for 
the CAP-TSDV calculations. This resulted in 
calculated static tip deflections of -1.33 in. at M= 0.628 
and +1.35 inch at M = 0.888. The Reynolds numbers 
for these two Mach numbers were 2.17 and 1.14 
million respectively based on the root chord of 2.0 feet 
for these test points in air. 

The calculated and experimental flutter boundaries are 
presented in Fig. 6. The experimental data show only a 
modest dip. The point at M=0.888 is estimated to be 
near the bottom of the transonic dip. Both the inviscid 
CAP-TSD and Euler results agree very well at the 
lowest Mach number but become significantly low or 
conservative at the highest Mach numbers. Although 

not shown, the linear theory calculations using CAP- 

42 

TSD by Gibbons ~ agree very well with the measured 
flutter data over the Mach range tested. When viscous 
effects are included using CFL3D-Navier-Stokes code, 
the calculated flutter points are much nearer the 
experimental boundary, but are still somewhat 
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conservative. The CAP-TSDV results are nearly the 
same as the CFL3D-NS results. 


The above computational results are for small 
amplitude transient analyses. A large amplitude limit 
cycle oscillation was calculated for the conditions 
shown as the solid symbols in Fig. 6 using CAP- 
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(a) Flutter speed index. 



(b) Frequency. 

Fig. 6 Comparison between experimental 
and calculated flutter speed index and 
frequency for a business jet flutter model 
tested in air. 


TSDV. Fig. 7 shows two transient responses 
illustrating this behavior. The motions were excited 
from converged statically deformed conditions by 
multiplying the modal displacements and velocities by 
5.0 for Fig. 7a and 0.5 for Fig. 7b. The large 
displacement gives a displacement of about 7 inches, 
which decays to about 5-6 inches peak-to-peak. The 
smaller initial excitation results in a similar 
oscillation. This is consistent with the observed flutter 
point at this Mach number of 0.888. The boundary 
layer calculations indicated intermittently separating and 
reattaching flow in the outboard regions. Increasing 
Reynolds number in the calculations resulted in a 
growing oscillation that led to a larger limit cycle. 




(b) Amplitude growing to limit cycle oscillation. 

Fig. 7 Calculated limit cycle response for a 
business jet wing flutter model. M = 0.888, 
q = 79 psf, Re c = 1.14 million. 
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Comments on Viscous Effects 


The results for the two wings discussed above have 
indicated strong viscous effects for some, but not all, of 
the cases. The 445.6 wings tested in air showed large 
calculated viscous effects on the backside of the flutter 
dip at supersonic speeds. The business jet was 
indicated to have large calculated viscous effects on the 
front side of the dip. The calculated limit cycle 
oscillations changed to a more divergent nature to a 
large limit cycle at a higher Reynolds number. 
Another case of viscous effects on flutter is indicated by 
the effect of transition strips on the flutter boundary of 
the Benchmark Supercritical Wing tested on the Pitch 
and Plunge Apparatus (PAPA) in ref. 43. The flutter 
boundary was raised about ten percent in dynamic 
pressure at low Mach numbers by adding a transition 
strip to the upper and lower surface of the model. 

In terms of anecdotal evidence, during flutter tests of a 
transport model in the TDT, the approach to flutter 
monitored using one over peak amplitude of wing 
response. This parameter was approaching zero 
systematically as dynamic pressure was increased 
indicating that a flutter point was near. However, as 
the expected flutter point was neared, tufts showed that 
the flow separated near the tip, and the oscillations 
became more random and buffeting-like. In this case 
the separation, which may have been related to tip 
aeroelastic twist, appeared to quench flutter. 

The role of viscous effects has been shown by a few 
isolated examples largely by inference from 

calculations. The generality of such effects is not yet 
clear. Although, it currently appears that the inclusion 
of viscous effects in the transonic range is essential, at 
least at model scale, the extrapolation of this trend to 
flight Reynolds Number is somewhat less certain. 
Current practice at the TDT is to use transition strips 
to assure a fixed transition location and a turbulent 
boundary layer. Further effort may be needed to study 
the best way to approximate viscous effects for flight 
Reynolds numbers at model scale. 


TRANSONIC SHOCK-BOUNDARY LAYER 
OSCILLATIONS 


Shock-boundary layer oscillations have been 
experimentally investigated for several airfoils. These 
oscillations are made up of separating and reattaching 
boundary layers coupled with moving shocks oftimes 
referred to as Shock-Induced-Oscillations or SIO. One 
notable example is the oscillations over the 18 percent 


circular-arc airfoil at transonic conditions and zero angle 
44 45-47 

of attack. Several investigators have been 

successful in calculating this flow. Included in the 
successful applications is the boundary layer coupling 
method used in CAP-TSDV? A comparison of 
measured and calculated frequencies using CAP- 
TSDV 37 is shown in Fig. 8. 



Fig. 8 Comparison of calculated and 
experimental SIO reduced frequencies for the 
18% circular arc airfoil. 
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Fig. 9 Comparison of calculated and 
experimental buffet onset boundaries for the 
NACA 0012 airfoil. 
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Another example is for the NACA 0012 airfoil at angle 
48 

of attack Oscillations on this airfoil occur at 

transonic speeds and moderate angles of attack, and 
mostly on the upper surface. The comparison of the 
experimental boundary with CAP-TSDV calculations is 
presented in Fig. 9. There have been other successes 
calculating this complex flow but some higher 
level codes do not capture these phenomena. Other 

novel nonlinear features of the SIO flows calculated 

37 

with CAP-TSDV are discussed by Edwards. 

These types of oscillations can also occur on 
51-52 

supercritical airfoils. The implications of such 

oscillations on three-dimensional flexible wings are not 

33 

clear. One brief test indicated that for a 3-D wing, 
the shock-boundary layer oscillations led to a buffeting 
condition and a low amplitude limit cycle oscillation in 
a higher vibration mode. 

CA CHALLENGES 

There are some data available that may assist in the 
evaluation of CA capabilities. A few are briefly 
discussed and involve various aerodynamic phenomena. 
In addition there are some difficult configurations that 
may need to be analyzed in the future and because of 
their difficulty may pose challenges to the state of the 
CA art. 

Rectangular Wings of Varying Thickness 

53 

A series of rectangular wings were tested to 
investigate the effect of thickness on flutter in the 
transonic speed range. The airfoil section was a simple 
circular arc and the series of wings varied in thickness 
from a flat plate to ten percent thick. Fig. 10 shows the 
flutter speed index for two of these wings, one with 4 
percent thickness and one with 8 percent thickness. 
Several investigators (ref. 54-57) have calculated 
specific cases, but the complete trend has not been 
demonstrated, possibly because it is a very large 
computational task. It has been found that the modal 
data must be calculated for the wings using plate 
elements. The flutter boundaries are also complicated 
by a moderate proximity of the divergence boundary to 
the flutter boundary. 

Benchmark Models Data 

The zero angle of attack flutter boundary for the 
Benchmark Model with the NACA 0012 airfoil is 



Fig. 10 Flutter boundaries for two 
rectangular wings with circular-arc airfoils 
of different thickness, ref. 53. 



Mach Number 


Fig. 11 Measured flutter boundaries for 
NACA 0012/PAPA model at zero angle of 
attack, ref. 58. 

shown in Fig. 11. The simple, rigid rectangular 

58 

wing is mounted on the Pitch and Plunge Apparatus 
(PAPA). The boundary to the left is the classical 
flutter boundary, which in this case rises with Mach 
number. However, near M=0.9, there is a sharp notch- 
like boundary that involves mostly plunging motion. 
At this condition, there are strong shocks near 0.75 
chord. Fig. 12 shows the boundary at M = 0.78 versus 
angle of attack. Above 5 degrees, the rapid decrease in 
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the flutter boundary is a type of shock-induced stall 
flutter with the flow separating and reattaching during 
the cycle of motion. These cases are for very high 
mass ratio and some unpublished calculations have 
shown that the classical boundary is sensitive to 
aerodynamic damping of the pitching mode. Two rows 
of pressure measurements are available. It should be 
noted that the stall flutter occurs near the shock- 
boundary layer oscillation region for the NACA 0012 
airfoil. 



1 

Amplitude 


o a = 0° 



Fig. 13 Dynamic response measurements for 
Aeroelastic Research Wing-2, q = 260-340 
psf, ref. 59. 

resorting to experiment. Another particularly notable 
case 60 is that of the B-2. 


Fig. 12 Measured flutter boundary variation 
with angle of attack for NACA 0012/PAPA 
model at M = 0.78, ref. 58. 

Aeroelastic Research Wing - 2 

During the wind tunnel test of the ARW-2 wing, large 

dynamic oscillations were encountered at Mach 

59 

numbers well above design conditions. Again it 
involved separated flows, and did not appear to be 
classical flutter-like or traditional buffeting. It is 
thought to involve movement of upper and lower 
surface shocks on the outer panel. The plot of 

1/Amplitude is given in Fig. 13. This is an 
aerodynamically complex situation on a configuration 
that is of modest configurational complexity. 

Limit Cycle Oscillations 

There are several aircraft that encounter limit-cycle- 
oscillations during flight. One semiempirical method 
that has had some success is described in ref. 2. The 
challenge is to calculate some of these cases using 
Computational Aeroelasticity techniques without 


Unusual Configurations 

There are many ongoing efforts to design more efficient 
aircraft through innovative configurations, such as the 
blended wing body, joined wings, and twin fuselage 
transports. An example of one such study is given in 
ref. 61. The geometric complexity of many of these 
configurations will be the formidable CA challenge. 
Although this paper has focused on the aerodynamic 
effects, many of the unusual configurations will have 
structural dynamic challenges as well. Integrated 
aerodynamic controls may become an integral part of 
the design and will also have to be accurately treated. 

CONCLUDING REMARKS 

The motivation for Computational Aeroelasticity and 
the elements of one type of the analysis or simulation 
process have been briefly reviewed. The need for 
streamlining and improving the overall process to 
reduce elapsed time and improve overall accuracy is 
evident. Further effort is needed to establish the 
credibility of the methodology, obtain experience, and 
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to build in experience for simplifying the method for 
future users. 

A variety of results for two flutter cases have been 
discussed. Reasonable agreement is shown for the 
codes when viscous effects are included. There is a 
compelling need for a broad range of flutter test cases 
for further comparisons as the current test cases used for 
evaluation of CA are few and limited in scope. Some 
existing data sets should be revisited with the current 
generation of CFD/CA codes and further data sets are 
needed for calibration efforts. 

Existing data sets that may offer CA challenges have 
been presented. These are the effect of thickness on 
flutter for rectangular wings with circular-arc airfoils; 
the classical flutter, shock-induced stall flutter, and 
plunge instability of the Benchmark Models wings on 
the pitch and plunge apparatus; and the buffet-flutter 
boundary for the ARW-2 wind tunnel model. Looking 
to the future, there will be significant challenges in 
calculating aeroelastic results for the many unusual 
configurations that are being considered. 
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